00 

o 

O 
(N 

00 



X 
Oh' 



ATLAS reach for Quarkonium production and polarization measurements 

Erez Etzion * and Jonatan Ginzburg ^ ^ on behalf of the ATLAS CoUaboration 

'^Tel Aviv University, Raymond and Beverly Sackler school of Physics and Astronomy. 
Tel Aviv 69978, Israel 

The ATLAS detector at CERN's LHC is preparing to take data from the first proton-proton coUisions expected 
in the next few months. We report on the analysis of simulated data samples for production of heavy Quarkonium 
states J/?/) and T, corresponding to an integrated luminosity of 10 pb~^ with center of mass energy of 14 TeV 
expected at the early ATLAS data. We review various aspects of prompt Quarkonium production at LHC: the 
accessible ranges in transverse momentum and pseudorapidity, spin alignment of vector states, separation of 
color octet and color singlet production mechanism and feasibility of observing radiative decays Xc and Xb decays. 
Strategies of various measurements are outlined and methods of separating promptly produced J ftp and T mesons 
from various backgrounds are discussed. 
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1. Introduction and motivations 

The LHC is expected to produce a large num- 
ber of quarkonium states such as J/ip and T al- 
ready in low luminosity runs during the first few 
years of running. These events are important for 
many physics studies. Their sizeable branching 
fraction into charged lepton pairs allows for easy 
separation of these events from the expected huge 
hadronic background. Being narrow resonances, 
the quarkonia arc perfectly suited for alignment 
and calibration of the ATLAS trigger and track- 
ing systems. On top of that, understanding the 
details of the prompt charmonia production is 
a challenging task and a good testbed for vari- 
ous QCD calculations, spanning both perturba- 
tive and non-perturbative regimes through the 
creation of heavy quarks in the hard process and 
their subsequent evolution into physical bound 
states. This note concentrates on the capabilities 
of the ATLAS detector to study the prompt char- 
monia production and polarization measurement 
at LHC. One approach which can be used to de- 
scribe the evolution of heavy quark antiquark pair 
into quarkonium bound state is the Color Evapo- 
ration Model (CEM) ^IJ. In CEM there is no cor- 
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relation of color and angular momentum quan- 
tum numbers between the initial QQ state and 
the final quarkonium state. It is assumed that 
soft gluon emission during the binding process 
alters the quantum numbers to arrive at the ap- 
propriate final state. Before the Tevatron experi- 
ments started to collect data it was assumed that 
the Color Singlet Model (CSM) [2] is the right 
explanation for quarkonium production mecha- 
nism. This model assumes that each quarko- 
nium state can only be produced by a QQ pair 
in the same color and angular momentum state 
as that quarkonium. Therefore, a J/Vj meson, 
for example, can only be derived from a cc pair 
created in a l^Si color-singlet state. The attrac- 
tive feature of the CSM is that it makes definite 
quantitative predictions for production cross sec- 
tions, because the binding probability is related 
to the quarkonium wavefunction from the poten- 
tial model. However, this approach predicted the 
J/iJj production rate an order of magnitude lower 
than that measured by CDF. The Nonrelativistic 
QCD (NRQCD) Colour Octet Model (COM) [5] 
was proposed to explain this discrepancy. In this 
efi'ective field theory QQ pairs produced with one 
set of quantum numbers can evolve into a quarko- 
nium state with different quantum numbers, by 
emitting low energy gluons. In contrast to the 
CSM, in which such transitions have zero prob- 
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ability. The good description of the Tevatron 
data by the COM model shown in Figure [1] is 
at least in part due to the tuning of some of it's 
parameters which were determined from the same 
data. However, the difference between the CDF 
measurement of 3/ip polarization dependence on 
transverse momentum (pr) and the theory pre- 
dictions (Figure [21) is motivating us to repeat the 
measurement in the LHC higher energy regime 
utilizing the advantage ATLAS will have in the 
number of events and the higher px reach. 



1 

0.8 
0.6 
0.4 
0.2 


-0.2 
-0.4 
-0.6 
-0.8 
-1 




10 



15 20 

Pj (GeV/c) 



25 



30 




5 10 15 20 

Pt (GeV) 



Figure 1. Differential cross-section of J/-0 pro- 
duction at CDF with theoretical predictions for 
color-singlet and color-octet model production [4] 



2. ATLAS di-muon trigger 

The LHC will collide two 7 TeV proton beams 
at a rate of 40 MHz, which together with pile-up 
will result in an interaction rate of up to around 
1 GHz. The ATLAS trigger acts to reduce this 
rate to around 100 Hz written to disks, whilst 
keeping only events 'of interest'. The trigger first 
level (LVLl) decision is based on coarse granular- 
ity of two sub-detector systems: the muon trig- 



Figure 2. Polarizations of production as a 
function of pr at CDF, (blue circles) , curves for 
limit cases of the kr factorization model, and a 
band for NRQCD predictions (from [5J). 



ger chambers and the calorimeters. The second 
level trigger (LVL2) receives from LVLl data re- 
stricted to limited Regions of Interest (Rol). For 
a LVLl muon, the LVL2 will use the informa- 
tion from the muon and inner detector tracking 
chambers to improve the muon momentum esti- 
mate, which allows a tighter selection based on 
this quantity. There are two specific types of trig- 
ger dedicated to quarkonium: one which requires 
two LVLl Rols corresponding to two muon can- 
didates with pt values above thresholds of 4 and 
6 GeV respectively, and the other requires a sin- 
gle LVLl Rol above a threshold of 4 GeV and 
searches for the second muon of opposite charge 
in a wide Rol at LVL2. We consider a com- 
plementary trigger on 3 ftp events triggering on 
a single higher px (10 GeV) muon and search- 
ing for a matched track at the off-line analysis 
level. Before incorporating trigger and recon- 
struction efficiencies, the predicted cross-sections 
for pp J/ip + were calculated for 

a number of px thresholds on the di-muon trig- 
ger (see Table [T] for details). Figures [3] illustrates 
the distribution of cross-sections across the values 
of the Pt of the harder and softer muon from the 
quarkonium decay without any muon cuts applied 
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and zero polarization. The lines overlaid on the 
plots represent the the muon trigger thresholds 
Pt > G GeV and > 4 GeV for the harder and 
the software muon respectively (denoted further 
by "6+4 GeV" or ^^6^l4), similarly "4+4 GeV" 
(/i4/i4) refers to the trigger threshold > 4 GeV 
applied on both muons and "10+0.5 GeV" {filO) 
refers to trigger pT threshold of 10 GeV applied 
on one of the muons only. In all cases a pseu- 
dorapidity, 77, of a muon lies within an interval 
hi < 2.5. 



the bulk of the production is in the region near 
the muon threshold of 5 and 4 GeV. This means 
that by lowering the threshold for the higher mo- 
mentum muon we can significantly increase the 
number of recorded T events. 

Estimated quarkonia cross-sections for the 
three scenarios are presented in Table [TJ The 
production will be dominated by J/ip and T(IS'). 
while V'' and T{2S), T(3S') are expected to give 
smaller contributions. 
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Figure 3. Density of J/tp production cross-section 
as a function of the hardest and softer muon px 
of muons coming from No cut was placed 

on the generated sample, but the overlaid lines 
represent the thresholds of observed events with 
trigger cuts applied. 



Table 1 

Predicted cross-sections for various prompt vec- 
tor quarkonium state production and decay into 
muons, with di-muon trigger thresholds /i4/Lt4 and 
^6^4 and the single muon trigger threshold /ilO 
(before trigger and reconstruction efficiencies). 
The last column shows the overlap between the 
di-muon and single muon samples. 



Quarkonium 




Cross-section, nb 
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28 


23 


23 
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T(15') 


48 
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T(2S) 


16 
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T(3S) 
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The charmonium production is composed of 
three main processes: direct singlet production, 
octet production and singlet production of x 
states. Each of these processes is characterized 
by different differential cross-sections, Fig- 
ure m illustrates the contributions of these three 
classes to the overall production rate for J/-0 once 
muon PT trigger cuts of 6 and 4 GeV are applied 
to the muons coming from the 3/ip. 

3. Charmonium reconstruction 



Even if the bulk of the 3/^ are produced un- 
der the trigger thresholds, nevertheless, due to 
the high cross-section the number of accepted J /ip 
will be substantially larger than at the Tevatron. 
In the T case, due to it's relatively larger mass 



In any event which passes the di-muon trigger, 
all the reconstructed muon candidates are com- 
bined into oppositely charged pairs, and each of 
these pairs is analyzed in turn. If the invariant 
mass of the two muons is above 1 GeV, we at- 
tempt to refit the tracks to a common vertex. The 
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Quarkonium production cross-sections 
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Figure 4. Monte Carlo diflFerential cross-section 
for J/ip production with J/tppr with the differen- 
tial cross-section contributions from direct color 
singlet (dashed line), cc states contribute to 3/tp 
through X decays (dotted line) and octet produc- 
tion (dashed dotted line) overlaid. 



fraction of surviving J/ip candidates before the 
refit is 96%. If a good vertex fit is achieved, the 
pair is accepted for further analysis. If the invari- 
ant mass of the refitted tracks is within 300 MeV 
for J /tit or 1 GeV for T (six times the expected 
average mass resolution) of the expected mass, 
the pair is considered as a charmonium candi- 
date. The reconstructed efficiency varies with pt 
and T]. When the J/ijj pT > 10 GeV we get a 
sharp rise in its acceptance. The rise in the T 
case is less sharp and it reaches a high accep- 
tance level around 30 GeV. Both channels reach 
a similar plateau at around acceptance of 85%. 
Due to the trigger requiring both muons to pass 
the trigger pT threshold, the angular separation 
between the two muons is not large for most of 
the accepted events. De scribing the opening an- 
gle by SR = -\- 5rp', where 5(j) and dr] are the 



differences of the azimuthal angle and the of the 
two muons from i/ip. Typical values for 5R are 
« 0.47 which means that requiring both muons to 
be above the px threshold forces them to fly very 
close one to each other. For that reason the J/V' 
angular acceptance follows closely the individual 
muon distribution and its dependence on mate- 
rial and detector effects. In the T case, due to its 
higher mass, it tends to be produced with higher 
Pt- For that reason the separation angle 5R is 
much wider, the muons do not go to the same 
area in the detector and therefore the efficiency 
dependence on pseudorapidity is much smoother 
than in the J/V' case. The main sources of low in- 
variant mass di-muons expected to dominate the 
background for prompt charmonium are: 

• Decays in flight of tt^ and - muons from 
this channel have a steeply falling momen- 
tum spectrum. Their contribution in the 
mass region of interest is expected to be 
very small where the requirement for ran- 
dom coincidences in the charmonium invari- 
ant mass range reduces it to below the 1% 
level. 

• Di-muon production via the Drell-Yan pro- 
cess - only small fraction of these events 
survive the di-muon trigger requirements, 
which makes this background essentially 
negligible. 

• Continuum of muon pairs from beauty (and 

charm) decays - contributions from charm 
decays have not been simulated, but de- 
spite having an estimated total rate twice 
as high as from beauty for a /i6/i4 trigger, 
as the Pt distributions of muons from charm 
quarks falls more steeply, charm sources are 
expected to contribute at a lower level than 
beauty. 

• Indirect i/tp production - this is the main 
source of expected background, its reduc- 
tion is discussed below. 

All the background sources apart from the Drell- 
Yan pairs contain muons which do not origi- 
nate from the interaction point, this is used to 
suppress their contamination by rejection of the 



ATLAS reach for Quarkonium measurements 

events containing a secondary vertex, if identified. 
Tfie radial displacement of the two-track vertex 
from the beamline is used to distinguish between 
prompt 3/ip and B-hadron decays having an ex- 
ponentially decaying pseudo-proper time distri- 
bution due to the lifetime of the B-hadrons that 
may decay into quarkonia. The pseudo-proper 
decay time, r, is defined as 

^ Pt{J/^) ' ^ ' 

where Mj/^ and pt{J/^) represent the 3/ip in- 
variant mass and transverse momentum, and L^y 
is the transverse decay length of the meson. The 
resolution in the pseudo-proper decay time is ex- 
pected to vary from 0.110 ps for the low pr char- 
monia down to 0.07 ps for the higher pT lighter 
angular dependence is very vague. As demon- 
strated in Figure [5] a cut on this quantity can ef- 
ficiently distinguish between prompt and indirect 
J/ip events. 

Using the ^10 trigger each reconstructed sin- 
gle muon candidate is combined with oppositely- 
charged tracks reconstructed in the same event 
within a cone of SR = 3.0. Any track, including 
those that were not identified as muons, are ex- 
amined. As in the di-muon analysis, we require 
that both the identified muon and the track are 
flagged as having come from the primary vertex. 
In addition, we impose a cut on the transverse im- 
pact parameter do, \do\ < 0.04 mm on the iden- 
tified muon and \do\ < 0.10 mm on the second 
track, in order to further suppress the number 
of background pairs from B-decays. We obtain 
a. J/tjj invariant mass resolution close to that in 
the di-muon sample. It's worth noting that the 
signal-to-background ratio around the 3/tp peak 
improves slightly with increasing transverse mo- 
mentum of J/V'- At higher pT the cos 9* accep- 
tance also becomes broader, which should help 
independent polarization measurements. 

4. Polarization studies 

The polarizations of quarkonia can be mea- 
sured using the angular distribution of the daugh- 
ter particles produced in the decay. 
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Figure 5. top - Pseudo-proper decay time, r, 
distribution for reconstructed prompt J/tp (cross- 
hatched, centred at zero) and the sum of prompt 
and indirect 3/iJj from B-decays (hatched, expo- 
nential distribution) . 

bottom - Sources of low invariant mass di-muons, 
reconstructed with a pdpA px trigger, with the 
requirement that both muons are identified as 
coming from a primary vertex and with a pseudo- 
proper decay time smaller than 0.2 ps. The white 
area represent the background that was rejected 
with the pseudo-proper decay time treatment. 
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— cx 1 + acos^ei*, (2) 

a cos 0* 

where cos 0* is the angle between the direction of 
the positive (by convention) muon from quarko- 
nium decay in the quarkonium rest frame and the 
direction of quarkonium itself in the laboratory 
frame. The polarization parameter a, defined as 
a — (<tt — 2(Tl)/{(Tt + 2(Tl), is equal to +1 for 
transversely polarized quarkonia production, (he- 
licity ±1). ctt and cfl are the transverse and lon- 
gitudinal cross sections. For longitudinal (helic- 
ity 0) polarization a is equal to —1. Unpolarized 
production consists of equal fractions of helicity 
states +1, and —1, and corresponds to a = 0. 

The previous Tevatron measurements were lim- 
ited to below around 20 GeV where the polariza- 
tion is best predicted and the theory most under- 
stood. At ATLAS we aim to measure the polar- 
ization of directly produced prompt quarkonium 
in the pT region up to ~ 50 GeV and beyond with 
extended coverage in cos 9*. This will allow for 
improved fidelity of efficiency, better discrimina- 
tion of longitudinal and transverse polarizations 
and therefore reduced systematics uncertainties. 

Two methods for polarization measurements 
have been examined. In the first one we used 
a templates chi-squared fit to the cos^9* distri- 
bution in six pt bins of Monte Carlo (MC)-data 
sample. A linear combination of three MC gen- 
erated templates histograms are fitted to the an- 
gular distribution of the data. The first is for a 
longitudinally polarized sample the second for a 
transversely polarized sample and the third is for 
background events sample. 

Feed-down from x states and B-decays dilute 
the prompt sample and lead to an effective depo- 
larization which is difficult to measure. The back- 
ground events sample was generated with zero 
polarization [a = 0) using the same detector ac- 
ceptance effects. The MC samples generated are 
large enough, and as a result the statistical fluctu- 
ations in the templates are negligible. The mea- 
sured values of a in six bins of pr are presented 
in Figure [H As one can see from this plot the 
average polarization is well concentrated around 
a — 0.5, which is the initial value of polarization 
that was set in MC-data sample. 
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Figure 6. The measured polarization parameter 
a and its error in six pT bin. The initial polariza- 
tion of a = 0.5 was set to the MC data sample. 
Statistics corresponding to integrated luminosity 
of about 5 pb~^. 



Evidently polarization measurement can signif- 
icantly suffer from low | cos 9* \ acceptance, and 
hence from difficulties in separating detector effi- 
ciency corrections from polarization state effects. 
The di-muon trigger requirement for both muons 
to be above a certain pT threshold (^6^4) reduces 
to minimum the acceptance at large values of 
I cos 9* I , where the difference between various po- 
larization states is the more pronounced (see for 
example the two polarization states on the lower 
plots of Figure [T]) . The acceptance of the single 
muon trigger sample (/.tlO) is very different. Here 
the efficiency is higher at large values of | cos 9* \ 
and drops in around zero. Figure [7] (on the top) 
demonstrates how the two samples complement 
each other mainly at the low px regions, while 
at high Pt the two triggers increasingly overlap, 
thus allowing for a cross-check of acceptance and 
efficiency corrections. In order to achieve this, the 
Pt distributions of both samples, /z6/x4 and /xlO 
events, were appropriately combined. The com- 
bined distribution dN'^'^'^ / d cos 9* , was corrected 
according to equation [3l 

diV™'' 1 dN'""^ 
^ . 

dcos9* A{pt, cos 9*) ■ El ■ 62 dcos9*' 
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Figure 7. Six figures at the top show kinematic 
acceptances of the /Lt6/x4 (solid red Hne) and filO 
(dashed blue lines) generator level cuts, calcu- 
lated with respect to the sample with no genera- 
tor level cuts on muon pT, in slices of J/f/j trans- 
verse momentum: left to right, top to bottom 
9-12 12-13, 13-15, 15-17, 17-21, above 21 GeV. 
Six bottom figures show combined and corrected 
(equation [3]) distributions in polarization angle, 
for longitudinally (dotted lines) and transversely 
(dashed hues) polarized J/ip sample in the same 
Pt slices with statistics corresponding to inte- 
grated luminosity of 10 pb 



Here ei stands for the trigger and reconstruc- 
tion efficiency, while £2 denotes the efficiency 
of background suppression cuts for each sam- 
ple, and A{pttC0s9*) is the kinematic accep- 
tance of the triggers selection. These cos 0* dis- 
tributions are fitted using equation [2l with a 
and normalization as free parameters for each pT 
slice. The results obtained when fitting unpolar- 
ized sample (a = 0) with statistics correspond- 
ing to 10 pb-i are: a = 0.156 ± 0.166, -0.006 ± 
0.032, 0.004±0.029, -0.003±0.037, -0.039±0.038 
and 0.019 ±0.057 corresponding to the pt bins as 
in Figure [T] The precision of the cross section de- 
rived from the normalization factor are ±4.35 in 
the first slice and decreasing from ±0.09 to ±0.04 
with increasing pT slices. Repeating the same 
study with T we get back numbers that are con- 
sistent with (a — 0) but with larger errors run- 
ning from ±0.17 to ±0.22. To further check the 
ability to measure the polarization we reweighted 
the raw MC distribution to emulate transversely 
polarized (a = +1) and longitudinally polarized 
(a = — I) J ftp samples. The same analysis was re- 
peated and the results can be seen on the bottom 
of Figure [7] 

5. Analysis of x production 

A sizeable fraction of prompt J /ijj and T are ex- 
pected to originate from radiative decays of heav- 
ier states, Xc and Xb- These states have even C 
parity and therefore have a strong coupling to 
the color-singlet two gluon state. About 30 to 
40 % of J/?/' in our signal will come from decays 
of Xc J/ip + 7- Unfortunately, the energies 
of these photons tend to be quite small. The 
ability of ATLAS to detect these photons and re- 
solve various x states is rather limited. Trying to 
evaluate that for each reconstructed charmonium 
candidate we calculated the invariant mass of the 
/i/Z7 system with all the photons found in that 
event. The fi^i-f system is considered to be a x 
candidate, if: the difference between the invariant 
masses of the fifij and /x/z systems lies between 
200 and 800 MeV, and the cosine of the opening 
angle between the J/V' and 7 is larger than 0.98. 

We fitted three guassians to the difference in 
invariant masses of the /i/i7 and /i/i measured 
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in those x candidates. The three mean posi- 
tions were fixed according to the 318, 412 and 
460 MeV corresponding to the xo, Xi and X2 
expected vahies. The fit parameters were the 
heights of the three Gaussian peaks h^Q, h^i, h^2, 
and the three parameters describing the smooth 
polynomial background of J/ip production from 
B-hadron decays which survive the pseudo-proper 
decay time requirement. The MC input am- 
plitudes of the peaks (15, 123 and 87, respec- 
tively) were reproduced reasonably well: h^o = 
15 ± iistat) ± lG{sys),h^i = 101 ± 4(stat) ± 
12(sys), h^2 = 103 ± A{stat) ± 9{sys) with strong 
negative correlation between the last two. The 
overall Xc reconstruction efficiency is estimated 
to be about 5%. 

6. Physics reach with early ATLAS data 

During the initial run of the LHC with the lu- 
minosity of 10^^ cm~^ s^"'^, one day of running 
corresponds to integrated luminosity of 1 pb~^. 
This can be translated to about di-muon trig- 
gered 15,000 J/V' fJ-y^ and 2,500 of T — > /x^. 
To that we can add the sample of 16,000 and 
2,000 J/V' and T triggered with one muon above 
10 GeV. In 10 days we expect to collect a sam- 
ple which is roughly equal to the MC statistics 
used in this note which may be enough for the 
first cross section and polarization pt dependence 
measurements. The precision of the i /ip polariza- 
tion can reach 0.02-0.06 depending on the level of 
polarization itself, as well as on our understand- 
ing the trigger and detector efficiencies, resolution 
and background estimations. 
By the end of one year, with the expected in- 
tegrated luminosity of 100 pb~^, the transverse 
momentum spectra are expected to reach about 
100 GeV and possibly beyond. With several mil- 
lion ^ /i/i decays, and better understanding 
of the detector, Xc JfpJ should become observ- 
able, while other measurements mentioned above 
will become increasingly precise. During the fu- 
ture high luminosity runs, we will have to increase 
the trigger threshold, and rescale the lower trigger 
threshold. Nevertheless the higher luminosity will 
further expand the range of reachable transverse 
momenta, and allow further tests of the produc- 



tion mechanisms and the polarization measure- 
ments, as well as make Xc reconstruction easier. 
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